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Abstract 



We present a scheme for determining if the quantum state of a small trapped 
atomic Bose-Einstein condensate is state with a well defined number of atoms, 
a Fock state, or a state with a broken U(l) gauge symmetry, a coherent state. 
The proposal is based on the observation of Ramsey fringes. The population 
difference measured by a Ramsey fringe experiment will exhibit collapse and 
revivals due to the mean- field interactions. The collapse and revival times 
depend on the relative strengths of the mean-field interactions for the two 
components and the initial quantum state of the condensate. 



I. INTRODUCTION 

Since the observation of Bose-Einstein condensation in trapped atomic gases in 1995 [|TJ, 
an unresolved issue in the theory of atomic Bose-Einstein condensates (BEC) is the qauntum 
state of the condensate H1|1|H. In the standard theory of BEC, which 
applies in the thermodynamic limit, the quantum state is one of well defined phase, (p. This 



state corresponds to a broken £7(1) gauge symmetry |H| [11]]. Because particle number, 



N, and phase obey the uncertainty relation ANAcf) > 1 |12[ , a state of well defined phase 
implies uncertainty in the particle number. 

Again, this result is based on assumptions which are valid only in the thermodynamic 
limit in which N — > oo. The second quantized Hamiltonian for a system of bosons is expressed 
in terms of the bosonic field operator, *&(r, £), and its adjoint which annihilate and create 
a particle at position r, respectively. The Hamiltonian is invariant under the global U(l) 
transformation, ^(r, t) — ► e ix ^/(r, t). The U (1) symmetry implies a conserved Noether charge 
which corresponds to the total number of particles. When a Bose-Einstein condensate is 
present, a single quantum state becomes macroscopically occupied and it is assumed that 
the field operator acquires a nonvanishing expectation value, 
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*(r,*))=V>(r,*), (1) 

with respect to the condensate state. Here, t) is the order parameter for the condensate. 
However, ip(r,t) is no longer invariant under the transformation, ip(r,t) — > e %x i\)(j, t) which 
implies that the U(l) symmetry is spontaneously broken in the condensate. Equation ([!]) 
may be rewritten as ($(r,t)^ = (<p(N - 1)| *(r,£) \<f(N)) where \ f(N)) and \<p(N - 1)) are 
"like" condensate states which differ by one particle Jl3| |14[]. In the thermodynamic limit, 
the difference between \tp(N)) and \(p(N — 1)) disappears, in which case the condensate is in 
a coherent state, if>(r,t) \ip(N)) = ip(r,t) \<p(N)) and ip(r,t) may be identified with the wave 
function for the quantum state in which Bose condensation has occurred (with the wave 
function normalization J d 3 r \i[)(t, t)\ 2 = N). However, it is not clear that Eq. ([!]) is still 
applicable when N is finite. Examples of BEC in condensed matter physics such as superfluid 
He may have N ~ 10 20 whereas trapped atomic gases typically have iV ~ 10 3 — 10 6 . 

In fact, there are two immediate objections to the use of a coherent state for finite particle 
number. First, at zero temperature, one expects that the true ground state of the condensate 
will be a number state (i.e. a Fock state) with no quantum fluctuations in the particle number 
even if we are ignorant of what that number is. Also, since the Hamiltonian for the system 
is U(l) symmetric, one must introduce a symmetry breaking field into the Hamiltonian in 
order to define the phase of the condensate [15|j. The symmetry breaking field vanishes in 



the thermodynamic limit. However, there is no physical interaction which corresponds to 
this symmetry breaking term and as such, it simply amounts to a mathematical trick [TO . 



The general definition of a Bose-Einstein condensate due to Penrose and Onsager [17 
is that the single particle density matrix, pi(r, r ; , t) = (&(r, t)^(r', t)\ does not vanish as 
for large separations, 

lim pi(r,r ,t) = $*(r, t)$(r', t). (2) 

|r— r'|^oo 

Although Eq. (fj) is consistent with Eq. ([[]), Eq. fl2|) simply requires the macroscopic 
occupation of a single quantum state and will therefore be true for a condensate that is in 
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a number state. Wright et al. have shown that, for small condensates (N ~ 10 3 ) described 
initially by a coherent state, the order parameter undergoes collapses and revivals such that 
ip(r,t) — > during the collapse but Eq. (fj) remains valid at all times since pi(r,r',t) is 
unaffected by the phase diffusion which causes the collapse and revivals ||. This implies 
that a coherent state description is inappropriate for small N since it is not an energy 
eigenstate of the system. Similar results were obtained in || where the depedence of the 
collapse and revival times on the dimensionality of the condensate and the trapping potential 
was studied. In contrast, Barnett et al. have argued that the best pure state description of a 
condensate is the coherent state since it is the most robust state with respect to interactions 
with the environment 0. In short, there appear to be no conclusive arguments for or against 
a coherent state description of atomic BEC's. 

Since a coherent state has a well defined phase, the appearance of interference fringes 
in the atomic density for two overlapping condensates would be an indication that the 
condensates were in coherent states (or some other superposition of number states such 
that one could ascribe a phase to the condensate). However, Javanainen and Yoo have 
shown that even if the two condensate states are initially in number states, there will be 
an observable interference pattern ||. This is because the destructive detection of atoms 
creates an uncertainty in the relative number of atoms in the two condensates since it 
not known from which condensate the detected atom came from. Consequently, with each 
atom detection, the relative phase between the two condensates becomes more precisely 
defined. Thus, any interference experiment based on destructive detection of atoms will not 
be able to distinguish between two condensates initially in number states or coherent states. 
Similar work has shown that the detection of spontaneously scattered photons between 
two condensates can establish a relative phase between the condensates even when the 
condensates are initially in number states (JTSJ . 



In this paper, we propose a method for distinguishing between a condensate that is in 
a number state and a coherent state. The method is based Ramsey's separated oscillatory 
field technique J19 in a two-component condensate such as 87 Rb [pip . Such an experiment 



has recently been performed by Hall et al. at JILA | 20f . Ramsey's technique, as applied to 



two-level atoms initially in their ground states, consists of applying two " |-pulses" generated 
by an external field of frequency u e which couple the ground state and excited state. These 
pulses are separated by a time T. The first pulse puts each of the atoms into a superposition 
of the ground and excited states with equal population. The relative phase between the two 
states then evolves as u Q T where hu is the energy difference between the two states. The 
second pulse creates a population difference between the two states which measures the 
relative phase accummulated by the atoms as compared to the phase accummulated by 
the external field during the period T. The population difference is then cos(5T) where 
5 = uj — u e . In a condensate, the population difference after the second pulse will be 
affected by two-body interactions which cause a phase diffusion of the relative phase of the 
two components in the interval T between the pulses. As such, the population difference 
will experience collapse and revivals as function of T. The collapse and revival times depend 
on the strength of the two-body interactions and the intial state of the condensate such that 
the collapse and revival times for a coherent state are different from a number state. 

Wright et al. predicted a similar effect for the interference fringe visibility of two spatially 
overlapping condensates ||. They showed that the revival time for the fringe visibility for 
condensates initially in coherent states was twice that of condensates in number states. How- 
ever, their result was based on the assumption that the intra-condensate interactions were 
the same and that inter-condensate interactions could be ignored. They also assumed that 
the coherence between the two number state condensates was established by measurement 
of the interference pattern in the same manner as described in [|J. 

The key advantage of the Ramsey fringe technique is that the collapse and revivals 
manifest themselves in the population difference between the two condensate components 
which is readily measured used absorptive or dispersive imaging of the condensate. Proposals 
to directly measure the order parameter, in order to detect collapse and revivals, such as |22| 



usually rely on the detection of scattered light from the condensate. Reference [^] involves 
two independent condensates that are in spatially separated potentials. Such an experiment 



would be technically difficult. 

The remainder of this paper is organized as follows: In section II we present the second 
quantized Hamiltonian for a two- component condensate and derive a two-mode model for 
the ground states of the two components. The two-mode Hamiltonian is then represented in 
terms of angular momentum operators by exploiting the equivalence between the algebra of 
two harmonic oscillators and the angular momentum algebra. In section III, we consider a 
condensate prepared in one of the modes with a state vector given by either a number state 
or a coherent state and study the time evolution of these states subject to two | pulses. 
In section IV and V, we discuss the collapse and revivals as well as relevant time scales for 
observing them. 



II. PHYSICAL MODEL 



A. Derivation of two-mode Hamiltonian 



We consider a collection of bosonic atoms that have internal states |1) and |2) with 
energies ftw /2 and —fkj /2, respectively. There is a spatially uniform time dependent 
radiation field with frequency u e which couples the two internal states with a Rabi frequency 
Q(t). The atom field detuning is denoted by 5 = uj — uj e . The atoms in states |1) and 
|2) are subject to isotropic harmonic trapping potentials Vi(r) =\m<jj(r 2 for % = {1,2}, 
respectively. Furthermore, the atoms interact via elastic two-body collisions through the 
interaction potentials Vij(r — r') = Uij5 (r — r') where Uij = 



— and a jo is the s-wave 



scattering length between atoms in states % and j. It is assumed that > corresponding 
to repulsive interactions. The Hamiltonian operator describing the system is given by, 



H — H atom + H co n 
H atom = [ d 3 rU{(r) 



ft 2 T r / \ 

-2^ V +yi(r)+ T 



*i(r) + *S(r) 



(3a) 
* 2 (r) 

(3b) 
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H coll = l - ! rfVjf/n^W^W^iW^iW+^^W^W^W^W 
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+2f/ 12 ^ t 1 (r)^ 2 (r)^ 1 (r)^ 2 (r)} . (3c) 

Here, H a tom is the single particle Hamiltonian and H co ii represents two-body interactions. 

The operators ^(r) and &j(r) are bosonic annihilation and creation operators for 
an atom in state i = {1,2} at position r which satisfy the commutation relations 
[^j(r),^t( r ')] = SijS(r - r') and [^(r), ^-(r')] = 0. The operators, *»(r), have been writ- 
ten in a field interaction representation which is rotating at the frequency of the external 
field, uo ei so that *i(r) = ^[ N \ r y^ an d vj> 2 ( r ) = #W( r )e-*" e */ 2 . Here, $f°(r) are the 
field operators in the normal representation. This explains the appearance of the detuning 
in Eq. (fj§). 

In the presence of the condensate, we assume that the field operators may be ap- 
proximated using a two-mode model such that ^i(r) = ai0i(r) + 5^i(r) and ^ 2 (r) = 
^202 ( r ) + 8^/ 2 (r) where the are the mode annihilation operators for the condensate 
modes which obey bosonic commutation relations a^, aj = Sij and [a«, aj] = 0. The 8^i(r) 
represent the field operator for the non-condensate modes and will be neglected since the 
number of atoms in these modes is assumed to be negligible compared to the condensate 



-£v 2 + ^(r) 



modes. For small condensates, such that Naij/aho,i < 1 where aho,i = J ^-is the har- 
monic oscillator length, <f>i(r) are given by the harmonic oscillator ground states of the trap, 
: (r) = ^p</>j(r) P3]| . Assuming a weak trap, aho,i ~ 10/im and ~ 5nm 
for 87 Rb, one has N < 2000. In the two-mode approximation the Hamiltonian becomes (with 

n=i), 

i i l - i - 

H — - (8 + 3ui) a\a x + - (-5 + 3u 2 ) a\a 2 + -Q*(t)a[a 2 + -Q{t)a\ai 

Jr \ (xi a i a i a i a i + X2d\a\a 2 a 2 + 2xi 2 a\aia\a 2 ^ ; (4) 



where 



n(t) = n(t) J d 3 r^(r)0!(r); (5a) 



Xi = U n J dVlMr)! 4 ; (5b) 
X2 = U 22 J rf 3 r|0 2 (r)| 4 ; (5c) 
Xi2 = U l2 J d 3 r|<Mr)| 2 |0 2 (r)| 2 . (5d) 

For Q(t) = 0, the eigenstates of Eq. (f|) are simply the number states \n\, n 2 ) F such that 
Ni \ni,n 2 ) F = rii \rii,n 2 ) F where iVj = ajaj. One may note that Eq. (||) with X12 = also 
describes tunnelling between two condensates in a double well potential . 

B. Angular momentum representation 

Equation (Q) may be expressed in a more convenient form by taking advantage of the 
mapping between the algebra for two independent harmonic oscillators and the algebra for 
angular momentum [^] . The mapping between the two algebras is achieved by making the 
following definitions 

J + = a\a 2 ; (6a) 
J_ = a\di\ (6b) 

Jz = 2 ( a ^ ai ~ <32 ' a2 ) ' ^ 

By noting that the x and y components of the angular momentum are given by the operators 
Jx = h(J+ + J-) and J y = ■k(J+ — J-), it follows that 

where N = Ni + N 2 is the total number operator which commutes with the two-mode 
Hamiltonian. Thus J 2 is a constant of motion with eigenvalues j(j + 1) and j = N/2. 
Consequently, for a state with definite N, Eq. (f|) has the form |J, 

H = AuJ z + X +J 2 Z + \to*{t)J+ + -£l(t)J- + - (xi + X2 + 2 X i2) f + \ (3^i + 3u 2 -xi- X2)j\ 

(8) 



where 

Aw = <J+|(w!-u; a ) + x-(2j-l); (9a) 

X+ = g (Xi + X2 - 2xi 2 ) ; (9b) 

X- = g (Xi - X2) ; (9c) 

In writing Eq. (|5J) we have made the replacement N — > 2j. This limits Eq. (^) to the 
subspace of states with the same total number of atoms. However, in the following section, 
we are only interested in calculating the expectation values of the operators J z and J± after 
some time t so that the replacement N — ► 2j does not affect any of our results even when 
we choose an initial state that is a superposition of number states. Therefore, we can drop 
the last two terms Eq. (H) so that our effective Hamiltonian is given by 

H = AujJ z + X+ J 2 Z + \w(t)J + + \mJ- (10) 

For fi(i) = 0, the eigenstates of Eq. (§) are simply the eigenstates of J z , \j, m) with m = 
—j, By noting that n\ = j + m and n 2 = j — m, it follows that \ j, m) = \j + m,j — m) F . 
In order to avoid confusion, we use the subscript F (F as in Fock) on the kets for the number 
state basis in order to distinguish them from the angular momentum kets. When Q(t) = 0, 
the condensate ground state is simply the lowest energy eigenstate of Eq. (g) for a fixed 
total number of atoms. This corresponds to the \j,m) state with m = integer ^— for 
integer ^— j < j and m = ±j otherwise. Here, integer() denotes the integer part of 
the number in parentheses. However, we assume that — < —j so that the ground state, 
\j, —j) = |0, N) F , is fully polarized. When the resonance condition 5 + | (u>i — dj 2 ) = is 
satisfied, \j, — j) will be the ground state when (i) \2 < X12 and x+ > or (ii) xi > X12 and 
X+<0. 

One might also consider a coherent state in the number state basis given by 

|0,a 2 ) = e"l-l 2 / 2 ^ | ,n 2 ) F = £ -^=«W I/, "iO . (") 



, v ™2 ! — n V^2 ! 

n 2 =0 v z n 2 =0 v z 
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as a variational wave-function that minimizes H with a 2 = yNe %tp . One may show us- 
ing Eq. (|) that 5E = (j, -j\ H \j, -j) - (0, a 2 \ H |0, a 2 ) = -X2N/2. Notice that X2 = 
U 22 / [(27r) 3/ X Di2 ] ~ U22/V where V is the volume of the trap. Therefore 5E ~ U 22 (N/V) 
which is an intensive quantity. Since (j, — j\ H \j, — j) is extensive, SE is negligible in the 
thermodynamic limit. Therefore, the coherent state represents a good variational wavefunc- 
tion for the the ground state energy in the thermodynamic limit. 

The primary advantage of the angular momentum representation is that the dynamics 
of the condensate can be understood in terms of a spin vector on a Bloch sphere. For strong 
external pulses of duration t p such that f* p dt \&(t)\ 3> |Au;| t p , \x+ \ t p , the time evolution 
operator, U = Te~ l f Hdt 1 is simply a rotation operator in spin-space 

R(9, 4>) = exp[— i6 (J x sin (ft — J y cos 0)], (12a) 

where 8 sin (ft = jj p dt Re (&(t)) and —8 cos (p = J* p dt Im . We have neglected the time 

ordering operator, T in Eq. ( |12a| ) . This is justified if Q(t) is a square pulse so that the Hamil- 
tonian commutes with itself at different times in the interval < t < t p . Equation ( |12a|) is 



a rotation in spin-space through an angle 8 about the rotation axis n = (sin0, — cos 0,0). 
When there is no external field present, the time evolution operator is simply 

U (t) = e -*(A W j,+ x+ J2)t (13) 
which is diagonal in the \j, m) basis. 

III. CONDENSATE DYNAMICS 

In this section, we consider the dynamics of the condensate subject to two external pulses 
separated by a time interval T. For t < we assume that there have been no pulses applied 
and that the condensate is in the ground state. We denote the two ground states, considered 
in the last section, at t = by 

\*n) = (14a) 
|*c)=|0,a 2 ); (14b) 
10 



These two states correspond to a spin vector that is pointing towards the south pole of the 
Bloch sphere. For the number state, \^n), the length of this vector in the — z direction is 
N/2. The coherent state, \^c) > also points in the —z direction but with an average length 
of N/2 and an uncertainty in the z-component of the length of AJ Z = \\fN ■ 

At time t = + , a pulse is applied that rotates the system about the y-axis through 
an angle of ~. This pulse is described by the rotation operator i2(~,7r). It has the effect of 
transferring half of the condensate population from state 2 into state 1 so that after the 
pulse (J z ) = with the spin vector now pointing in the +x direction. The condensate is 
then allowed to evolve freely for a time T according to Eq. (0). After this period of free 
evolution, a second |-pulse is applied, again given by the rotation -R(§, 7r). After the second 
pulse, the population difference, (^Ni — iV 2 y = 2 (J Z (T)) as a function of T is measured. 

For T = 0, the effect of the two pulses would just be a spin- flip so that the spin vector 
would now be pointing in the +z direction. Because the state of the system following the 
first pulse (for both \^n) an d l^c*)) is a superposition of \j, m) states for all m values, 
the free evolution due to Eq. ([□]) causes the spin vector to diffuse in the equatorial plane 
as the different \j, m) states get out of phase with each other. Because the m are discrete 
integers, the \ j, m) states can re-phase leading to a revival of the spin vector. This dephasing 
and rephasing of the spin vector manifests itself as collapse and revivals of the population 
difference following the second pulse. However, the two initial states, \^n) and l^c); have 
very differenct collapse and revival times owing to the fact that \^c) is a superposition of 
states with different j. 

The calculation of (J Z (T)) is straight forward and we outline the calculation for the two 
intial states in the following two subsections. 

A. Number state, |\l/jv) • 

The quantity we wish to calculate is 

(UT)) N = {*n\ R{\, 7r)t[/ (T)tit:(|, 7r)t J,J2(E, 7r)U (T)R(^, tt) . (15) 
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First we note that 

R(P, n) j J z R((3, 7r) = cos (3J Z - sin (3J X (16) 
so that -R(f, 7r)t J z R(~, 7r) = — Re{J + }. Consequently, Eq. (|T5|) reduces to 



(Jz(T)) N = - Re { <tf „| J?(|, 7r)t[/ (T)t J + [/ (T)i2(|, tt) |tf „>} = - Re{(J + (T)) N } (17) 



The matrix elements of R(^,n) are easily calculated for arbitrary j |2l)|| , so that the state 
of the system following the first pulse and the free evolution period is simply, 



m=-j 



(18) 

Finally, one obtains 

(J + (T))n = -je^cos 2 ^ 1 ( X+ T) ; (19a) 
so that the population difference, (^Ni — N 2 ^ = 2 (J Z (T)) N is then 

Ni - N 2 } = Ncos ((d+ | ( Wl - cj 2 ) + x_(iV - 1)) Tj cos^ 1 ( X+ T) . (20) 



B. Coherent State, \V C ) ■ 

The calculation of (J Z (T)) C is similar to that of ( J Z (T)) N ,the main difference being an 
average over a Poissonian distribution of number states. As before, we only need to calculate 
the expectation value of J + following the free evolution period, 

(J,C0> o = -Re{<¥ o |i2(|,^ (21) 
Following the first pulse and the free evolution period, the state of the system is 
U (T)R{~, tt) \Hf c ) = e-l-l 2 / 2 £ (22) 



n 2 =0 v z 



X 



J^. —y—y U-m)\U + rn)\ ^ H + ^ 
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Again, it should be emphasized that Eq. ( |22|) is only valid for calculating matrix elements 
for operators that are diagonal in j, which include all angular momentum operators. The 
evaluation of (J + (T)) C can be done in two steps, 



00 1 |2n 2 
|2 \ a 2 ' 



(J + (T)) C = e-l-l 2 £ ^V<W, (~je iA » T cos 2 ^ 1 ( X+ T)) , (23) 



Tl2=0 

where the term in parantheses is the same as Eq. ( |19aj ). The final step is just an averaging 
over a Poissonian distribution of the number of atoms (|ct 2 | 2 — N), 

(J+(T)) C = ~e< 6+ ^- u 2) ) T exp [N (-1 + cos ( X+ T) e^~ T )] . (24) 

The population difference following the second pulse is then, 

(Wi-iV^ = ATexp[Ar(-l + cos(x + T)cos(x-T))] (25) 
x cos f h$ + I ( Wl - w 2 )^ T + AT cos (x+T) sin(x-T)J . 

Equations ( p0|) and (|25| ) represent the central results of this paper. Note that although 
the number of atoms, N, appearing in the two equations has the same value, the meaning 
of N is different. In Eq. (pC|), N is the exact number of particles while for Eq. ( p5|) , AT is 
the average number of particles for a superposition of number states. 



IV. COLLAPSE AND REVIVALS 

From Eqs. ( f20|) and (|25|) one can see that the population difference involves a rapidly 
oscillating part and an envelope function that is responsible for the collapse and revival of 
the population difference. For simplicity we assume that the external field is on resonance 
so that 5+| {uj\ — uj 2 ) = and the population difference for the two cases simplify to 

(aTi-A^ =Ncos(x-(N-l)T)cos N ' 1 ( X+ T); (26) 
<^A>i - N 2 } = Nexp [N (-1 + cos (x+T) cos(x_T))] cos (N cos (x+T) sin(x-T)) ; (27) 

In the following sub-sections we consider several limiting cases. 

13 



A. x- = 



The simplest nontrivial case to consider is Xi = X2 X12 (this corresponds to 



Ref. f§ where X n = 0) so that ^JV X - N 2 J = iVcos 7 ^ 1 (x+T) and - N 2 
iVexp [N (— 1 + cos(x + T))] . This case is shown in Figure 1. The population difference 
quickly decays to zero for both cases as soon as cos(x+^) deviates significantly from 1. 
The collapse time may be estimated by making a Gaussian approximation for small times. 
One finds then for x+T <C 1 

(Ni -N 2 ^j ~ Ne- {T/ ™ )2 ; (28) 
Nj. ~N 2 ^ c ~ Ne- {T/Tc)2 ; (29) 



and the collapse times are tn = X+ 1 V^/ ffi ~ 1) an d T c — X+ 1 ^/^/^ which are indistin- 
guishable for N 1. The variance in J z following the first pulse is given by 



AJ Z = y/(^|it;(-,7r)J|it;(-,7r) |^) = ^/J/2 = VN/2 (30) 

for i = N,C. The collapse times can be expressed as r N = r c = 1/ [\f2x+AJ z ) which 
shows that the collapses are attributable to the dephasing of the different J z states due to 
the x+ Jz term in the Hamilitonian. 

The revival times are quite different for the two states. For the number state, the revivals 
occur whenever T/v = mr/x+ where n is an integer. However, when the number of atoms is 
even, iV — 1 is odd and the condensate will undergo anti-revivals when T = (2n + 1)tt/x+ 
so that (^Ni — iV 2 ^ = —N at these times. (Note that a similar affect was described in 
Ref. in which the fringe visibility of the interference pattern could undergo a revival 
with a 7r phase shift when the number of atoms that had been detected in order to establish 
an interference pattern was even). In contrast, the coherent state undergoes revivals at 
the times Tq = 2mr/x+ which is twice the revival time of the number state. In addition, 
^JVi — iV 2 ^ > regardless of N. Therefore, there appear to be two key differences that 
distinguish a number state from a coherent state for X- — 0: (i) The occurrence of a negative 
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population difference when N is even and (ii) revival times that are half the revival times 
of a coherent state. 

Since the revivals are determined by the x+^1 i n Eq. (0), there is a simple explanation 
for the factor of two difference in the revival times. The state U (T)R(^,n) \^n) involves 
a superposition of | j, m) in which the m values all differ by an integer. Consequently, the 
revivals occur when the relative phase between all of the \j, m) states is an integer multiple 
of 2n. This corresponds to the condition [(m + l) 2 — m 2 ] x+T = 2n'ir + for all m and 
where is a global phase factor that is independent of m and n' is an integer. By taking 
n! = mn and = x+T, one sees that the revivals occur at integer multiples of the time tc/x+- 
However, for the coherent state one has instead U (T)R(j, 7r) \^?c) which is a superposition 
of states with different j and m values so that in this case the values of m in the superposition 
need only differ by a half-integer. Therefore, for the coherent state the condition for the 
occurrence of a revival is [(m + 1/2) 2 — m 2 ] x+T = 2n'ir + 0. Again taking n! = mn but 
with = x+T I '4, the revivals occur at integer multiples of 2n/x+- 



B. x- + and X-N > 1 



In this case Eqs. (|26l - p7|) consist of a rapid oscillations modulated by a slowly varying 



envelope function which gives rise to the collapse and revivals. As such, we only consider 
the behavior of the envelope functions in this sub-section which is given by Jn{T) and fc(T) 
for the number state and coherent state, respectively. The envelope functions are given by 

f N (T) = cos^- 1 (x+T) , (31a) 
fc(T) = exp [N (—1 + cos ( X+ T) cos(x-T))] . (31b) 

The collapse and revival times for the number state are the same as what was found in the 
previous subsection. The only difference is that there are no antirevivals since Jn{T) = — 1 
simply corresponds to a n phase shift in the rapidly oscillating part of Eq. 

However, the behavior of the coherent state is quite different. One can see that for 
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short times, fc(T) ~ e~^ T / TC " >2 where the collapse time is Tc = \J w^T^ )' Therefore, 
increasing \- decreases the collapse time. This is attributable to the X-(2j — 1)J Z term in 
the Hamiltonian which causes the states with different j but the same m to get out of phase 
with each other in a time ~ l/(x_A7V) = l/(x-VN)- The reduction in the collapse time is 
illustrated in Figure 2 for the case X-/X+ — 2. 

Revivals occur when cos (x+T) cos(x_T) = 1 which can only be satisfied if X-/X+ — p/q 
where p and q are integers. When this condition is satisfied, the revivals occur at times 
Tc = n7i/x+ where n is a positive integer and (p/q)n is an integer such that if n is odd 
(even) then (p/q)n is also odd (even). When X-/x+ is irrational, there are no revivals and 
even for rational values of X-/x+i the period between revivals can differ significantly from 
the revival period for the number state, 7r/x+. For example, if X-/x+ = 1/4 then the first 
revival will occur at 8n/x+ for the coherent state. In Figure 3, the revivals are shown for the 
ratio x~/x+ = 1/3- One can see that the number state has revivals at all integer multiples 
of n/x+ while the first revival for the coherent state occurs at 3n/x+- 



V. DISCUSSION 

In the previous section it was shown that the population difference in a Ramsey fringe 
experiment can exhibit collapses and revivals with times that can be very different for a 
number state and a coherent state. For a number state, we will, in general, be ignorant of 
the exact number of atoms so that we can only say that the state \0,N) occurs with some 
probability p(N). Consequently, an ensemble average over many different experimental runs 
will lead to statistical fluctuations in the number of atoms even though for the number 
state, the number of atoms in any given experiment would be precisely defined. Suppose 
p(N) obeys a Poisson distribution, will there be any difference now between the number 
states and the coherent state? There is in fact a difference since the quantum fluctuations of 
the coherent state are present in each experimental run whereas the statistical fluctuations 
of the number states only become manifest after averaging over many such runs. As an 
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example, consider the case of the anti- revivals for the number state, ^Aq — N 2 j = —N, 
which can occur when xi — X2 7^ Xi2- In this case, one might observe such anti-revivals in 
any single experiment with a number state but such anti-revivals would never be observed 
in any experiment for a coherent state since it is the averaging over the Poisson number 
distribution which prevents the anti-revivals. Even though the anti-revivals may be ob- 
served in each experimental run for the number states, the average over many experimental 



runs, ^2 N p(N) ^Aq — N 2 j = \Ni — N 2 j , will not exhibit the anti-revivals since they get 
averaged out. 

One may also consider initial states of the condensate given by |^) = J2 n c « \®-> u )f where 
c n is sharply peaked around n = N (such clS Cji OC e (™ Ar )/ 4cr ). This state satisfies Eq. (P 
for finite N. However, will exhibit qualitatively similar behavior to the coherent state M 
since the critical difference in the results of the Ramsey fringe experiment lie in the quantum 
fluctuations in the particle number of the coherent state. 

Finally, we estimate the order of magnitude of the collapse and revival times. The 
fundamental time scale which determines the collapses and revivals is H/x+ (where we now 
explicitly include that factors of H). If the trapping potentials for the two components are 
the same, uq = u 2 = u, then, 0i(r) = 4> 2 {r) = <f>(r) and J \(f)(r)\ 4 d 3 r = 1/ [(27r) 3//2 a| ] where 
dho = \Jh/muj. If we assume a relatively weak trap, a^o = 5jjm and if we take U%j = ~ 



A-!rh 2 a 



with a = 5nm as an estimate for 87 Rb, one finds Xi/h ~ (^7^) / [(2 7r ) 3 ^ a L] 



0.023s _1 . Therefore the collapse times will be on the order of {Xi/fy~ ~ Is for 

= 1000 and the revival times will be on the order of tt {Xi/fc)" 1 ~ 100s. Note that since 
the scattering lengths are nearly equal for the two hyperfine states of the 87 Rb condensate 
P 1|| , x± ~ for U\ = uj 2 . However, this may be overcome by either manipulating one of the 
scattering lengths using a Feshbach resonance PB| or by changing the trapping potentials 
so that (f>i (r) 7^ 4>2{r). 

As mentioned, a Ramsey fringe experiment has been performed at JILA ||20|| . However, 
this experiment was performed with a relatively large condensate (N = 5 x 10 5 ) such that 
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Naij/ciho,i 3> 1 [23|. The theoretical analysis of the experiment |2(| |^7| was based on 
the Gross-Pitaevskii equation which is a mean-field equation for the order parameter and, 
as such, already assumes a state of broken symmetry in the condensate. In addition, the 
duration of this experiment (i.e. the period between pulses, T) was too short too observe 
the phase diffusion due to the quantum dynamics of the condensate (see footnote [23] in 
Ref. |E0|). 



VI. CONCLUSION 

In this paper we have shown that Ramsey's separated oscillatory field technique applied 
to a small atomic Bose-Einstein condensate exhibits collapse and revivals in the population 
difference between the two internal states of the condensate. The collapse and revival times 
depend on the strength of the two-body interactions and the initial state of the condensate 
so that one may potentially distinguish between a condensate state that is a number state 
or a coherent state. Since absorptive and dispersive imaging of atomic BEC measure the 
density of atoms in the condensate, the Ramsey fringe experiment proposed here should 
be easier to perform than an experiment which tries to directly observe the collapse and 
revivals in the order parameter. 
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Figure 1. Plot of the normalized population difference, (^Ni — N 2 j /N as a function of 
T for the case X- — an d N = 1000 for the number state, (a), and the coherent state, 
(b). 

Figure 2. Plot of the normalized population difference, (^Ni — N 2 ^ /N for x+T <C 1, 
X-/x+ — 2, and N = 1000 illustrating the rapid oscillations and the envelope function 
which leads to the collapse. The collapse time for the coherent state is significantly 
smaller than the number state. 

Figure 3. Plot of the normalized population difference, (^Ni — iV^ /N as a function of 
T for X-/x+ — 1/3 and N = 1000 which shows that revivals at multiples of tt/x+- ( a ) 
is the number state and (b) is the coherent state. 
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